First-principles study on the intermediate compounds of LiBH4 
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We report the results of the first-principles calculation on the intermediate compounds of LiBH4. 
The stability of LiBsHg and Li2B„H„(n — 5—12) has been examined with the ultrasoft pseudopoten- 
tial method based on the density functional theory. Theoretical prediction has suggested that mon- 
oclinic Li2Bi2Hi2 is the most stable among the candidate materials. We propose the following hy- 
driding/dehydriding process of LiBH4 via this intermediate compound : LiBH4 <^ j^Li2Bi2Hi2 -I- | 
LiH -f y|H2 «-+LiH -f B -)-|H2. The hydrogen content and enthalpy of the first reaction are estimated 
to be 10 mass% and 56 kj/mol H2, respectively, and those of the second reaction are 4 mass% and 
125 kJ/mol H2. They are in good agreement with experimental results of the thermal desorption 
spectra of LiBH4. Our calculation has predicted that the bending modes for the F-phonon fre- 
quencies of monoclinic Li2Bi2Hi2 are lower than that of LiBH4, while stretching modes are higher. 
These results are very useful for the experimental search and identification of possible intermediate 
compounds. 



I. INTRODUCTION 

Hydrogen is the most promising environmentally clean 
energy carrier to replace fossil fuels. The use of hydrogen- 
based energy in practical applications such as fuel cell ve- 
hicles, however, requires the development of safe and effi- 
cient hydrogen storage technology. Complex hydrides in- 
cluding the light metal lithium (Li) have sufficient gravi- 
metric hydrogen storage capacity, and many researches 
and development about the lithium complex hydride such 
as LiBH4^'^'^ and LiNH2'*"'''*' have been done recently. 
Particularly, LiBH4 can desorb about 14 mass% of hy- 
drogen by the following thermal decomposition : 



LiBH4 



LiH + B + -H2. 



(1) 



However, the experimental value of enthalpy of this reac- 
tion is 69 kJ/mol H2^ and LiBH4 is too stable to release 
hydrogen at ambient condition. High temperature and 
high pressure are needed for rehydriding reaction, so its 
reversibility becomes a problem for practical use. 

Ziittel et alX investigated the hydrogen desorption 
from LiBH4 in details, and reported that the thermal des- 
orption spectra of LiBH4 mixed with Si02-powder exhib- 
ited three hydrogen desorption peaks. These peaks were 
observed around 500 K, 550 K and 600 K, and the prod- 
ucts corresponded to "LiBHa.g" , "LiBHs" and "LiBH2" 
, respectively. The compounds in quotes were nominal 
compositions estimated from amount of desorbed hydro- 
gen. This is a strong evidence for the existence of the 
intermediate compounds , and the hydrogen desorption 
reaction takes place in at least two steps not a single step 
like in Eq.(^. Although low temperature release of hy- 
drogen and improvement of reversibility can be expected 



by use of the intermediate compound, little attention has 
been paid yet. 

In this study, the stability of the intermediate com- 
pounds of LiBH4 has been investigated theoretically and 
we clarify a hydriding/dehydriding process of LiBH4. 
The thermal desorption experiment observed a structural 
transition around 380 K and the melting of the compound 
at 550 K. The stability of the solid phases at absolute 
zero temperature is mainly discussed here. In Sec.m we 
describe the details of the computational method and 
the possible intermediate compounds of LiBH4. Sec- 
tion mil reports the calculated results on the stability 
of some candidate compounds. And then, we discuss 
the hydriding/dehydriding reaction via the intermediate 
compounds. Furthermore, the electronic structures and 
T-phonon frequencies on the most stable compound are 
studied. We also investigate the stability of the borane 
complex anions, and then compare with the solid com- 
pounds. 



II. COMPUTATIONAL DETAILS 

First-principles calculations have been performed by 
the ultrasoft pseudopotential method^ based on the den- 
sity functional theory J. The generalized gradient approx- 
imation (GGA) formulck^ is applied to the exchange- 
correlation energy. 

The interaction between the ion cores and electrons 
is described by the ultrasoft pseudopotential^, and the 
norm-conservation constraint-P^ is imposed on Li for 
the calculation efficiency improvement. The scalar- 
relativistic all-electron calculations are first carried out 
on the free atoms(ions). We chose 2s and 2p states for 
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both Li and B pseudopotential as the reference states 
with the cutoff radii of 2 a.u. (Li) and 1.5 a.u. (B), re- 
spectively. A single projector function is used for each 
angular momentum component. The 3d state is treated 
as the local part of pseudopotential. The hydrogen pseu- 
dopotential is constructed from Is, 2s and 2p states with 
the cutoff radii of 1.1 a.u.(s) and 1.2 a.u. (p). We 
use double projector functions for the s component and 
a single projector function for the p component. For 
all pseudopotentials, the pseudo-wave functions and the 
pseudo-charge-augmentation functions are optimized by 
a method similar to that proposed by Rappe et ali^ Also, 
the partial core correctioni^ is taken into account for Li 
and B pseudopotentials. 

In the solid-state calculations, the pseudo-wave func- 
tions are expanded by plane waves with a cutoff energy 
equal to 15 hartrees. The cutoff energy for the charge 
density and potential is set to be 120 hartrees. The inte- 
gral over the Brillouin zone is approximated by the sum- 
mation on the k-grid elements of which the edge lengths 
closer to the target value of 0.15A~-'^ as possible. We 
confirmed that these calculation conditions gave a good 
convergence of energy within 0.002 eV/atom. The pre- 
conditioned conjugate-gradient technique is employed to 
minimize the Kohn-Sham energy functional. A proce- 
dure based on the iterative diagonalization schemei^ and 
the Broyden charge mixing methodic is adopted in this 
study. Optimization of crystal structures is performed till 
the atomic forces and the macroscopic stresses become 
less than 5 x 10~^ hartree/bohr and 0.1 GPa, respec- 
tively. During the structural optimization process, the 
partial occupation numbers near the Fermi level are de- 
termined by the Fermi-Dirac distribution function with 
ksT = 3 X 10^'^ hartrees. The Helmholtz free-energy 
functional, including the entropy term, is minimized 
instead of the Kohn-Sham energy functional. Then, the 
improved tetrahedron method^^ is used in order to mini- 
mize the Kohn-Sham energy functional in the optimized 
structure. The dynamical matrix is calculated by the 
force-constant method'''^ to obtain the F phonon frequen- 
cies. The atomic displacement is set to be 0.02 A. The 
further details of calculation are described in Refs. 0, 
and references therein. 

As the candidates of the intermediate compounds, the 
existing alkali metal-B-H materials are used. For exam- 
ple, it is well known that boron (B) and hydrogen (H) 
form inorganic compounds called 'borane' and the com- 
pounds of CsBgHg^ KsBgHsjia and K2Bi2Hi22a are 
reported. The space group of crystal structure for the 
compound CsBsHg is Ama2 (No. 40) and cation Cs"^ and 
anion [BaHg]^ are arranged same as the NaCl-type struc- 
ture with an orthorhombic distortion. The crystal struc- 
tures of K2BeH6 and K2B12H12 are best described as an 
anti-CaF2-type arrangement with K"*" cation in the cen- 
ter of a tetrahedron formed by [B„H„]^~ dianions. We 
first assumed that LiBsHg, Li2B6H6 and Li2Bi2Hi2 com- 
pounds, which had the same crystal structures as exist- 
ing CsBsHg, K2B6H6 and K2B12H12, respectively. The 



stability of these candidates was evaluated using the first- 
principles calculations. Since the existence of a series of 
closo-type dianions [B„H„]^"(n = 5 — 12) is also well 
known, our calculation have been expanded to Li2B„H„ 
(n = 5 — 12) compounds. Although Li atom and [B„H„] 
cluster are arranged in the anti-CaF2-type structure and 
the unit cell parameters are supposed to be a face- 
centered-cubic symmetry (a = b = c,a = P = S = 90° 
) as starting points for the structural optimization, the 
output relaxed compounds have different symmetry de- 
pending on the structure of [B„H„] clusters. The al- 
kali metal salt of monoclinic K2Bi2(OH)i2i2i with closo- 
[Bi2(OH)i2]^^ dianions has been reported, and mono- 
clinic Li2Bi2Hi2 with similar crystal structure is also ex- 
amined. 

III. RESULTS AND DISCUSSIONS 
A. Stability of candidate compounds : LiBa,Hy 

Table U shows the calculated results on the struc- 
tural parameters and the cohesive energies of LiBsHg 
and Li2B„H„ (n = 5 — 12). We denote the cubic 
Li2Bi2Hi2 based on K2B12H12 with type-1 and mon- 
ochnic Li2Bi2Hi2 based on K2Bi2(OH)i2 with type- 
2, respectively. Compared with the cohesive energy 
of two type of Li2Bi2Hi2, the value of monoclinic 
Li2Bi2Hi2(type-2) is larger than that of type-1. There- 
fore, the type-2 Li2Bi2Hi2 is easy to form. After the 
next paragraph, the only result concerning monoclinic 
Li2Bi2Hi2 (type-2) is shown. 

The enthalpy of formation for hydriding reactions in- 
cluding LiBH4 from LiH (NaCl-type), a-B (rhombohe- 
dral) , and H2 molecule are given in Table (nj where the 
zero-point energy corrections are not taken into consider- 
ation. They are provided using calculated cohesive ener- 
gies of 2.3609 eV/atom for LiH, 6.2013 eV/atom for a-B, 
2.2689 eV/atom for H2 molecule, and 3.1501 eV/atom 
for LiBH4 (orthorhombic Pnma symmetry), respectively. 
The enthalpies of formation for Li2B„H„ (n = 10, 11, 12) 
are more negative than that for LiBH4. Therefore these 
compounds have a great potential for generating in hy- 
driding reactions from LiH and a-B as the intermediate 
phase of LiBH4. 

In order to understand the stability of each compound 
intuitively, we represent the enthalpy of formation cor- 
responding to the following reaction with hydrogen of 6 
mole in Fig. ^ : 

(l-^)[LiH+B]+(5H2 ^-^UB^Ky+^—^^^^—^UE, 

X X 

(2) 

where 5 — {y — l)/(2x -|- y — 1) shows the mole fraction 
of H2 and a; > 1. From Fig. ^ monoclinic Li2Bi2Hi2 
is formed as the intermediate compound of the hydrid- 
ing/dehydriding reaction of LiBH4, because the actual 
reaction goes along the lowest state of enthalpy. There- 
fore, following hydriding/dehydriding process is proposed 
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TABLE I: Structural parameters and cohesive energies {Ecoh) of the candidate compounds : LiB^Hy. It is denoted that the 
cubic Li2Bi2Hi2 based on K2B12H12 with type-1 and monochnic Li2Bi2Hi2 based on K2Bi2(OH)i2 with type-2, respectively. 



Compound 


Space group 


Unit cell parameters 


Atom 


Site 


Atomic position 

x/a y/b 


z j c 


Ecoh 

(cV/atom) 


LiB3 Hg 


Ama2 


a ^ 9.188 A 


Li 


46 


0.0595 


0.6321 


0.25 


3.3770 




{No. 40) 


b — 8.813 A 


Bl 


46 


0.1372 


0.1335 


0.25 








r> — ^ 7fi'^ A 
C — . / Do 


B2 


8c 


0.3090 


0.1331 


0945 










HI 


Ah 


0.0690 


0.2507 


25 










H2 


46 


0.0684 


0.1666 


25 










H3 


8c 


0.1386 


0.5149 


0.5190 










H4 


8c 


0.1823 


0.1333 


0063 










H5 


8c 


0.3619 


0.2508 


0.0186 




1j12 05 -n.5 


RSm 


„ — c cnq A 


Lil 


3a 








0.2755 


3.7910 






(- — 1 fi 76*^ A 


Li2 


3a 








0.7513 






(hexa.goii£il axes) 




Bl 










0.0876 










B2 


3a 








0.9295 










B3 


96 


0.1070 


0.8930 


0.0084 










HI 


3a 








0.1595 










H2 


3a 








0.8586 










H3 


96 


0.2270 


0.7730 


—0.0078 




1j12 Do -n-D 


If to< f t 


n — 7 QfiS A 


Li 


8c 


0.25 


0.25 


0.25 


4.0429 




(iVo.225) 




B 


24e 


0.1553 
















H 


24e 


0.3054 










Li B H 

2 7 7 


/2 


a = 5.691 A 


Lil 


26 





0.1501 


0.5 


4.0850 




(iVo.5) 


X, _ Q QOfJ A 
u — ;?.oou 


Li2 


26 


0.5 


0.2317 











^ — ^ fi^f^ A 

c — 0.0^0 >^ 


Bl 


2a 





0.1564 


Q 








(3 = 77 39° 


B2 


4c 


-0.1871 


0.0529 


Q 1063 










B3 


4c 


0.1160 


-0.1087 


0659 










B4 


4c 


-0.1335 


0.0087 


0.1903 










HI 


2a 





0.2792 













H2 


4c 


-0.3463 


0.0805 


-0.2099 










H3 


4c 


0.2215 


-0.2030 


0.1343 










H4 


4c 


-0.2700 


0.0109 


0.3879 




1j12 jj8 -H-S 


/42m 


^ — K ^79 & 

Ui — O'OI^ -rt 


Li 


4d 





0.5 


25 


4.0920 




TVrt 1911 


C — XU.OOl A 


Bl 


8i 


0.1032 


0.1032 


Q \24:2 








B2 


Si 


0.1631 


0.1631 


0304 










HI 


Si 


0.2031 


0.2031 


— 0.2112 










H2 


Si 


0.3114 


0.3114 


0.0580 




T j1 n rTn 


R3 


a — 7,044 A 


Lil 


3a 








0.2482 


4.1326 




( Mr, 1 


^ — 1 c rifi9 A 
C — i-O.UOZ 


Li2 


3a 








8070 






(hexagonal axes) 




Bl 


96 


0.8087 


0.9031 


u.youo 










B2 


96 


0.8118 


0.9047 


0.0694 










B3 


96 


0.8661 


0.7287 


0.0096 










HI 


96 


0.6730 


0.8346 


0.8923 










H2 


96 


0.6756 


0.8357 


0.1281 










H3 


96 


0.7699 


0.5322 


0.0141 




1 ji ri 1 c\ ti 1 n 

J-J-i^ J-' J. U -1- u 


7422 


a = 6.196 A 


Li 


id 





0.5 


0.25 


4.2360 




{N0.97) 


c = 10.356 A 


Bl 


4e 








0.1827 










B2 


16A: 


0.1942 


-0.0809 


0.0739 










HI 


4e 








0.2984 










H2 


16/c 


0.3647 


-0.1582 


0.1070 




1j12 JJii rill 


12 


« — 7 ns^ A 


Lil 


26 





-0.1091 


0.5 


4.2397 




yri 0.0) 




T 


OK 


0.5 


-0.1083 


u 








c = 7.U0U A 


Bl 


2a 





-0.1264 











/3 = 70.52 


B2 


4c 


0.1117 


0.1190 


— 0.1120 










DO 


4C 


0.0835 


-0.0162 


— 0.2427 










B4 


4c 


0.2425 


-0.0162 


n nszin 

— . u ^ u 










Pit; 


4C 


0.1293 


0.0825 


1290 










B6 


4c 


0.1448 


-0.0795 


1446 










HI 


2a 





-0.2355 













H2 


4c 


0.2016 


0.2052 


—0.2020 










H3 


4c 


0.4160 


-0.0284 


—0.1910 










H4 


4c 


0.1902 


-0.0283 


— 0.4165 










H5 


4c 


0.2133 


0.1471 


0.2127 










H6 


4c 


0.2305 


-0.1441 


0.2300 




i-jl2 -t> 12 ^ 12 


Frn3 


a = 10.083 A 


Li 


8c 


0.25 


0.25 


0.25 


4.2997 




{N0.202) 




B 


ASh 





0.1448 


0.0881 










H 


48h 





0.2490 


0.1469 




J-jl2 012 J^12 


P2\j 71 


n — 7 A 

ti ( .000 rt. 


Li 


4e 


0.6747 


0.6250 


u . 00^0 


4 3461 




( ATn 1/11 


— Q ^^P. A 
— ?J.OOD A 


Bl 


4e 


0.6770 


0.5392 


1632 








c — 6.768 A 


B2 


4e 


0.4447 


0.5629 


0.2239 








/3 = 92.26° 


B3 


4e 


0.5425 


0.6760 


0.0442 










B4 


4e 


0.6915 


0.5689 


0.9022 










B5 


4e 


0.6856 


0.3922 


0.9959 










B6 


4e 


0.5329 


0.3887 


0.1933 










HI 


4e 


0.8033 


0.5642 


0.2786 










H2 


4e 


0.4008 


0.6050 


0.3840 










H3 


4e 


0.5707 


0.7983 


0.0742 










H4 


4e 


0.8252 


0.6151 


0.8271 










H5 


4e 


0.8231 


0.3231 


0.9941 










H6 


4e 


0.5558 


0.3097 


0.3310 
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TABLE 11: The enthalpies of formation for the hydriding re- 
action of various compounds LiB^^H^ including LiBH4 from 
LiH, a-B, and H2 molecule, where the zero-point energy cor- 
rections are not taken into consideration. 



Hydriding reaction 


Enthalpy of formation 
(kJ/mol H2) 


LiH 4- 3B + 7/2H2 


LiBsHg 


-36 


2LiH + bB + 3/2H2 


Li2B5H5 


113 


2LiH 6B 2 H2 


^ Li2B6H6 


-42 


2LiH + 7B + 5/2H2 


— > Li2B7H7 


-45 


2LiH 4- 8B -)- 3 H2 


— + Li2B8H8 


-32 


2LiH -f 9B + 7/2H2 


Li2B9H9 


-42 


2LiH -f lOB 4 H2 


— > Li2BioHio 


-87 


2LiH + IIB + 9/2H2 


—> Li2BiiHii 


-79 


2LiH + 12B + 5 H2 


— > Li2Bi2Hi2 


-125 


LiH + B + 3/2H2 


^ LiBH4 


-75 



-3 

Q. 
CD 
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FIG. 1: Enthalpy of formation relative to (1 - 5) [LiH + B] 
and molecular H2 of 5 mole (reaction Eq.l|5J in text), as a 
function of the mole fraction of H2 , S. 



LiBH4 ^ — Li2Bi2Hi2 + -LiH+ — H2, (3) 
12 6 12 ' 

^ LiH + B + ^H2. (4) 

The enthalpy without zero-point energy effects and 
hydrogen content of the first reaction (Eq.||3I)) are 
56 kJ/mol H2 and 10 mass% and those of the second re- 
action (Eq.(0J) are 125 kJ/mol H2 and 4 mass%, respec- 
tively. These results agree with the experimental ones, 
which reported that 9 mass% of hydrogen are liberated 
in the hydrogen desorption peak in thermal desorption 
spectra of LiBH4 mixed with Si02-powder in Ref. 1. The 
enthalpy value of 56 kJ/mol H2 for the first reaction is 
lower than computational one of 75 kJ /mol H2 for direct 
dehydriding reaction of LiBH4 (Eq.(^)i^ So the hydrogen 
desorption and absorption can occur at low temperature 
and pressure comparatively by using this intermediate 
compound. 



We also calculated other possible crystal struc- 
tures such as Cu2BioHio-typ922, and [Li(thp)3]2[BiiHii]- 
typei^ Moreover, we performed the calculation for the 
F-phonon frequencies of the candidate materials. If there 
was a soft-mode, we lowered crystal structure symme- 
try with moving the atoms along the direction of the 
soft-mode eigenvectors. The maximum gain of cohesive 
energy was 0.03 eV/atom for Li2B6Hg, however, no com- 
pound or crystal structure more stable than monoclinic 
Li2Bi2Hi2 was found. The present calculation result on 
the energy of monoclinic Li2Bi2Hi2 provides the upper 
limit value for the thermodynamic stability. In con- 
clusion, the existence of the intermediate compound of 
LiBH4 was predicted theoretically. 

Recently Kang et al.^'^ have reported that LiBH and 
LiB are the intermediate phases of LiBH4, and propose 
the dehydriding reaction of LiBH4 through LiBH as fol- 
lows : 



LiBH4 



LiBH+ -H2. 
2 



(5) 



We also performed the first-principles calculation of the 
orthorhombic phase of Pnma LiBH. The reaction en- 
thalpy for Eq.|Hl) is 1.30eV per LiBH4 formula unit 
(84kJ/molH2), which is in good agreement with reported 
value of 1.28eV per LiBH4 formula unit in Refl23. On the 
other hand, the reaction enthalpy of Eq.Q is 0.63eV per 
LiBH4 formula unit. Therefore, the dehydriding reaction 
via the intermediate compound Li2Bi2Hi2 is energeti- 
cally more preferable one. 



B. Intermediate compound : Li2Bi2Hi2 

Here we describe the fundamental properties, such as 
the crystal structure, the electronic structure, and the F 
phonon frequencies, on Li2Bi2Hi2 which is expected as 
intermediate compound of LiBH4. 

Optimized crystal structure model of monoclinic 
Li2Bi2Hi2 is shown in Fig. [3 The bond lengths be- 
tween B and H atom of Li2Bi2Hi2 are 1.20 - 1.21 A, 
and they are very close to those for LiBH4 (1.23 — 1.24 
A) reported in Ref. |3- As for the B-B bond lengths 
of intra-icosahedron, the values of 1.779 — 1.811 A for 
Li2Bi2Hi2 are comparable to the experimental ones for 
the a-rhombohedral boron (a-B) of 1.751 — 1.806 A,— 
too. Since boron crystal has the icosahedral B12 cluster 
as a common structural component, the decomposition of 
[Bi2Hi2]^~ anion into B12 cluster and hydrogen molecule 
is easy to understand. 

Figure 13 shows the total and partial density of states 
(DOS) for Li2Bi2Hi2, which denotes that it has the en- 
ergy gap of 5.60 eV. Since there is little contribution of 
Li orbital for occupied states, Li2Bi2Hi2 consists of Li+ 
and [Bi2Hi2]^~ ions. The orbitals of B and H hybridize 
each other and the feature of occupied states in DOS is 
analogous with the distribution of the computed energies 
of the bonding molecular orbitals in [Bi2Hi2]^~i2& 
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FIG. 2: (Color online) Crystal structure model of monoclinic 
Li2Bi2Hi2(type-2). Large, middle, and small spheres denote 
Li, B, and H atoms, respectively. 



The F-phonon mode frequencies of monoclinic 
Li2Bi2Hi2 have been calculated so that the vibrational 
properties can be compared with experiments easily and 
directly. The phonon density of states is shown in Fig. 01 
It is divided into three regions in the same case as LiBH4. 
The first region is less than 300 cm^^ where the displace- 
ments of Li atoms are dominant. The second region is 
between 450 cm~^ and 1100 cm~^ where the B-H bond 
of [Bi2Hi2]^~ dianions vibrates with changing in angle 
between them (bending modes), and the third region 
is between 2400 cm~^ and 2600 cm~^ where the inter- 
atomic distance of B-H bond is changing along bond axis 
(stretching modes). The frequencies of bending modes 
are lower than LiBH4 described in Ref. 0, while stretch- 
ing modes are higher. 

The investigation using in situ Raman spectroscopy is 
effective for the confirmation of the short-range order or 
bonding of LiBH4.'^^ We examine the atomistic vibrations 
of LiBIl4 during heating by in situ Raman spectroscopy, 
and the identification of spectra modes originating from 



Li2Bi2Hi2 is now in progress 
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C. Stability of complex anions : [B„H„]^ 

Finally, we consider the relation between the stability 
of the compounds Li2B„H„ and those of complex an- 
ions [B„H„]^~. The energies of the isolated [B„H„]^~ 
are obtained using a face-centered-cubic supercell with 
a = 20 A. The single F point is used for the k-point 
sampling. The energies for the charged systems are com- 
puted by adding uniform background charges and im- 
proved with Makov and Payne correction^^ for the inter- 
action between periodic image charges. 

Figure |S1 shows the comparison of the formation ener- 
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FIG. 3: (Color online) The total and partial densities of 
states (DOS) for Li2Bi2Hi2. The energy is measured in eV 
relative to the top of valence states. 
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FIG. 4: Phonon density of states for monoclinic Li2Bi2Hi2. 
The contribution of the TO P-phonon modes indicated by ver- 
tical bars is only taken into account and the Gaussian broad- 
ening with a width of 30 cm~^ is used. 



reactions: 



n - 2 

2LiH + nB-\ — H2 ^ Li2B„H„ 



and 



?iB + -H2 + 2e" 



[B„H„] 



(6) 



(7) 



where the energies are normalized by the number of B-H 
pairs for comparison purposes. We can find a fairly good 
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FIG. 5: Comparison of the formation energies, ET'"" and 



-icomplex 



corresponding to the reactions Eq.Q and Eq.Q in 
text, respectively, where the energies are normalized by the 
number of B-H pairs. 



correlation between both energies. This is probably due 
to the fact that the electrostatic interaction between Li^ 
and [B„H„]^~ is not sensitive to n and a large part of this 
energy cancels out with that of LiH in Eq. Among 
the c/oso-type dianions considered here, [Bi2Hi2]^~ is the 
most stable one. 

LiBH4 desorbs hydrogen at temperatures above the 
melting point. The latent heat of fusion has been re- 
ported^- to be 0.078 eV/formula-unit for LiBH4 and 
similar values are expected for Li2B„H„, which is con- 
siderably smaller than the energy difference between 
[Bi2Hi2]^~ and other cZose- type dianions. The stability 
of the [Bi2Hi2]^~ anion supports our main conclusion in 
this study, that is, the intermediate phase suggested by 
the experiment is Li2Bi2Hi2. 



ing hydriding/dehydriding process of LiBH4 is proposed 
: LiBH4 ^ l/12Li2Bi2Hi2 + 5/6LiH -fl3/12H2 ^ LiH 
+ B + 3/2H2. The hydrogen content of the first and 
the second reaction are 10 mass% and 4 mass%, respec- 
tively, which agree well with the thermal desorption spec- 
tra (TDS) experiment on LiBH4.^ The heat of formation 
without zero point energy corrections for the first reac- 
tion, which is estimated from the solid state LiBH4 with 
Pnma symmetry (not liquid LiBH4), is 56 kJ/molH2. 
This value is lower than that for the direct reaction 
(LiBH4 ^ LiH + B +IH2). Therefore, low tempera- 
ture release of hydrogen can be expected by use of this 
intermediate compound. 

We have calculated the electronic structure and the F- 
phonon frequencies of monoclinic Li2Bi2Hi2. This com- 
pound has the energy gap of 5.60 eV and consists of Li+ 
and [Bi2Hi2]^^ ions. From the phonon density of states, 
it is predicted that its bending modes have lower fre- 
quencies than that of LiBH4,^ while stretching modes 
are higher. The identification of the experimental Ra- 
man spectra modes originating from Li2Bi2Hi2 is now in 
progress. 

The stability of closo borane complex anions [B„H„]^~ 
was also examined. We found a fairly good correlation 
between the formation energies of the solid phases and 
the isolated dianions. This result supports the validity 
of the intermediate compound indicated in the TDS ex- 
periment being Li2Bi2Hi2. 

There are various kinds of borane, such as nido-type 
[B„H„4_4] and arac/ino-type [B„H„+6], in addition to 
c/oso- type borane [B„H„_|_2]^ The nido and arachno bo- 
rane are derived from closo borane by removing one and 
two vertices, respectively. These generate the salts in the 
ionized state as well as closo-type dianions. These alkali 
metal salts are also the candidates of the intermediate 
compound of LiBH4. We will study the details of the 
hydriding/dehydriding process for LiBH4 including the 
stability of these materials in future. 



IV. SUMMARY 

We have examined the stability of LiBsHg and 
Li2B„H„(n = 5 — 12) , which are the possible interme- 
diate compounds of LiBH4, by the first-principles cal- 
culation. Our computational results for enthalpy of the 
hydriding reactions provide that monoclinic Li2Bi2Hi2 is 
the most stable one among the candidates. The foUow- 
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